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We investigate the properties of domain wall fermions on a set of quenched configurations at non-zero temper- 
ature. In particular, we compute the low lying eigenvalues of the DWF operator and study their relation with 
topology, level crossings and chiral symmetry breaking. We also measure the screening correlators and discuss 
the dependence on the extent of the extra dimension and the quark mass. 



1. Introduction 

An interesting aspect of the Domain Wall 
Fermion (DWF) formulation is that it provides a 
consistent framework for studying the chiral limit 
Tdq — > 0. In particular, it is expected to satisfy 
the Atiyah-Singer index theorem. This is a dis- 
tinct advantage over the Wilson and staggered 
formulations which are well known to suffer from 
shifts of the low-lying eigenvalues in the real and 
imaginary directions respectively 0]. However, 
the DWF formulation accomplishes this feat at 
the expense of introducing a 5 th dimension and, 
strictly speaking, all the nice properties are only 
recovered in the limit N5 — > oo (overlap formal- 
ism). A careful (non-perturbative) study of the 
convergence is therefore needed before practical 
simulations can be carried out. In this study, we 
initiate such a program in the context of QCD at 
non-zero temperature (thereby continuing our ef- 
fort to understand the link between topology and 
the chiral phase transition 0]). We have used 
for our computations a set of quenched configu- 
rations on a 16 3 x 8 lattice. The (3 values [and 
number of configurations] studied are: 6.2[170], 
6.1[170], 6.0[100], 5.9[100] and 5.8[100]. These 
configurations were used earlier to compute the 
screening correlators and the low- lying eigenvalue 
spectrum with a staggered fermionic action [Q. 
Here, we first study the level crossings of the 
hermitian Wilson operator (in a manner similar 
to H), then move on to compute the spectrum 
of the DWF operator (also properly hermitized) . 
The DWF calculations are carried out (so far) 



at N5 = 4, 6, 8 and 10. The link between topol- 
ogy, the number of level crossings and the number 
of DWF eigenvalues which vanish (exponentially) 
wiht N§ can then be studied. As was noted ear- 
lier for staggered fermions 0], the knowledge of 
the low-lying eigenvectors is also quite valuable 
since in the high temperature phase, their contri- 
butions dominate the disconnected mesonic cor- 
relators. We confirm this again for DWF and use 
it to improve the convergence of our conjugate 
gradient inverter and to lower the error bars on 
the disconnected correlators. 
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Figure 1. Level crossing of the hermitian Wil- 
son operator on a configuration with topological 
charge 1. 
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Figure 2. Eigenvalues of the hermitian DWF 
operator vs. A5 on a configuration with Qtop = 1. 



2. Results 

In this preliminary report, we focus our at- 
tention on the results obtained at (3 — 6.2 (the 
highest temperature we have studied) . The com- 
plete results, including the trends observed as one 
lowers the temperature towards the chiral phase 
transition will be reported elsewhere ||. Also 
because of lack of space we illustrate the vari- 
ous computations on a single configuration. We 
have selected a configuration which has topolog- 
ical charge 1. We then obtain one level crossing 
in the spectrum of the hermitian Wilson opera- 
tor (Fig. 1). It occurs at /1 « 0.78 with a slope 
of 0.74. On our sample of configurations, most 
crossings occured between fi = 0.75 and fi — 0.80. 
One notable exception was the single configura- 
tion of topological charge 2 for which the second 
crossing occured at fi ~ 0.95 . There is a relation- 
ship between the presence of these crossings and 
the existence of zero modes of the hermitian op- 
erator H = 75 i?D (where D is the DWF operator 
and R is the reflection operator in the 5 th dimen- 
sion). For example, a crossing at /1 with slope 
±1, would imply an exact zero- mode of H at a 
domain wall mass M = 1 + fi. Around this value 



and for slopes of magnitude less than 1 , there will 
be near zero-modes. In this study we have fixed 
M = 1.7 (similarly to what has been done in zero 
temperature applications [Q). On the same con- 
figuration as before, We then find that the lowest 
eigenvalue of H decreases exponentially with N$ , 
whereas the second lowest is essentially indepen- 
dent of N 5 (Fig. 2). The lowest eigenvalue is of 
order 0.0016 at N5 — 10 and would be 2 order of 
magnitudes smaller if A5 was extended to 20. 
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3.2 : conf. § 24 lowest eig. 
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Figure 3. 4-dimensional slices of the "Left" and 
"Right" components of the lowest eigenmode 



Our variational method for computing the 
eigenvalues also gives us the eigenvectors. Fig. 
3 represents the positive and negative chirality 
components measured on the lowest mode as a 
function of x§ . The negative chirality component 
is peaked on the left wall (as should be for a con- 
figuration with Qtop — —1) and decreases expo- 
nentially in the interior. The right chirality com- 
ponent also reaches its maximum on the left wall 
but with a smaller magnitude. For modes with 
no net chirality, one would obtain a symmetrical 
graph with each component peaked on their re- 
spective wall. (This is what is observed for the 
"first excited state" on our sample configuration 
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and also for the lowest mode on configurations 
with Qt p — 0). On a few "abnormal" config- 
urations with very late crossings ||, we found 
eigenmodes with a net chirality but an eigenvalue 
which remains large at A5 = 10. A detailed dis- 
cussion of this (small) artifact is left for || . 

Finally, as an introduction to our computa- 
tion of disconnected mesonic correlators ||, we 
present our results for the measurement of m 
Tr^$S as a function of m at N 5 — 10. The con- 
tinuum answer should be equal to the topological 
charge of the configuration ( - 1 in our case ) . At 
low values of the quark mass, fig. 4 shows that the 
measurements dip towards 0. This is as expected 
for a situation where the lowest eigenvalue is not 
quite yet at N5 = 10. It is also important to 
note that a precise measurement of this opera- 
tor was made possible by "separating out" the 
contribution from the lowest eigenmodes in the 
conjugate gradient inversion. Good results could 
then be obtained by using only 10 noise vectors 
on each wall. Without this "trick" , even increas- 
ing the number of noise vectors by a factor of 10, 
still leaves us with rather large error bars (see Fig. 
5). 
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Figure 4. Chirality versus quark mass (using our 
knowledge of the low-lying modes) 



3. Conclusions 

We have confirmed (non-perturbatively) that 
DWF provide a practical and systematically im- 
provable way of studying the topological prop- 
erties of QCD (at least in the high temperature 
phase). The "rate of convergence" at (3 — 6.2 
can be read off from Fig. 2 (and similar results 
at lower (3 will be presented elsewhere). In ad- 
dition, we have shown how the knowledge of the 
low-lying eigenmodes can be used to greatly im- 
prove the quality of fermionic measurements on 
configurations with non-trivial topology. We also 
expect that enhancements of this technique will 
play an important role in the design of dynamical 
fermion algorithms at very low quark masses. 
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Figure 5. Chirality versus quark mass (noisy 
estimator only) 
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